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' A reduction of the mass of the 7/(958) meson may indicate the restoration of the Ua(1) 

symmetry in a hot and dense hadronic matter, corresponding to the return of the 9th, 
"prodigal" Goldstone boson. We report on an analysis of a combined PHENIX and STAR 
^ | data set on the intercept parameter of the two-pion Bose-Einstein correlation functions, as 

I . measuremed in ^snn = 200 GeV Au+Au collisions at RHIC. To describe this combined 

PHENIX and STAR dataset, an in-medium r?' mass reduction of at least 200 MeV is 
needed, at the 99.9 % confidence level in a broad model class of resonance multiplicities. 
Energy, system size and centrality dependence of the observed effect is also discussed. 
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Although the quark model exhibits a U(3) chiral symmetry in the limit of massless up, down 
and strange quarks, and in principle 9 massless Goldstone modes are expected to appear when 
this symmetry is broken, only 8 light pscudoscalar mesons are observed experimentally. This 
puzzling mystery is resolved by the Adlcr-Bcll-Jackiw Ua{1) anomaly: instantons tunneling 
between topologically different QCD vacuum states explicitely break the Ua(1) part of the U(3) 
symmetry. Thus the 9th Goldstone boson is expected to be massive, and is associated with the 
rj' meson, which has a mass of 958 MeV, approximately twice that of the other pseudoscalar 
mesons. 

In high energy heavy ion collisions at RHIC, a hot and dense medium is created. Recent 
measurements of the direct photon spectrum in = 200 GeV Au+Au collisions indicate pQ , 

that the initial temperature in these reactions is at least 300 MeV, while hadrons as we know 
them may not exist above the Hagedorn temperature of Tjj ~ 170 MeV [2]. Thus the matter 
created in heavy ion collisions at RHIC is hot enough to be a quark-gluon plasma [T] . Detailed 
analysis of the properties of this matter indicate that it flows like a perfect fluid [3] , and scaling 
properties of the elliptic flow indicate scaling with the number of constituent quarks @], hence 
this matter is sometimes referred to as a strongly interacting Quark-Gluon Plasma (sQGP) [5], 
or, in more direct terms, a perfect fluid of quarks [3]. 

After this perfect fluid of quarks rehadronizes, a hot and dense hadronic matter may be 
created, where the Ua{^) symmetry of the strong interactions may temporarily be restored 
[7J IB]- Recent lattice QCD calculations indicate that such chirally symmetric but hadronic 
matter may exist below the critical temperature for quark deconfinement [3] • In such a medium, 
the mass of the 77' (958) mesons may be reduced to its quark model value of about 500 MeV, 
corresponding to the return of the "prodigal" 9th Goldstone boson [7] . In this presentation we 
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summarize the the results on an indirect observation of such an in-medium 77' mass modification 
based on a detailed analysis of PHENIX and STAR charged pion Bose-Einstein correlation 
(BEC) data (TUJOJl- These results have been published recently in Refs. [T3l [14] and detailed 
in Ref. [15], 

The abundance of the 7/ mesons with re- 
duced mass may be increased at low pt, by 
more than a factor of 10. One should empha- 
size that the 77' (and 77) mesons almost always 
decay after the surrounding hadronic matter 
has frozen out, due to their small annihilation 
and scattering cross sections, and their decay 
times that are much longer than the charac- 
teristic 5-10 fm/c decoupling times of the fire- 
ball created in high energy heavy ion collisions. 
Therefore one cannot expect a direct observa- 
tion of the mass shift of the 77' (or 77) mesons: 
all detection possibilities of their in-medium 
mass modification have to rely on their en- 
hanced production. 

An enhancement of low transverse momen- 
tum 77' mesons contributes to an enhanced pro- 
duction of soft charged pions mainly through 

the 7/ -> 1] + 7T+ + 7T~ -> (lT + + 7T° + 7T~) + 

n + + 7r~ decay chain and also through other, 
less prominent channels. As the 77' decays far away from the fireball, the enhanced production 
of pions in the corresponding halo region will reduce the strength of the Bose-Einstein corre- 
lation between soft charged pions. The transverse mass (ror = y m 2 + p\) dependence of the 
extrapolated intercept parameter A* of the charged pion Bose-Einstein correlations was shown 
to be an observable that is sensitive to such an enhanced 77' multiplicity, as pointed our first in 
Ref. [16] and discussed in Ref. [IT] . 

The predicted drop of A*(mx) data at low transverse masses has been observed both by 
PHENIX QJl HU and STAR [12 [18] at RHIC, as it is indicated on Fig. [fl Note that this 
reduction is not present in the S+Pb data set at CERN SPS at y r SMN = 19.4 GeV energy [19]. 

2 Modeling and analysis method 

Our main analysis tool was a Monte-Carlo simulation of the transverse mass dependence of 
the long lived resonance multiplicities including the possibility of an enhanced rj' production at 
low transverse momentum, due to a partial in- medium Z7a(1) restoration and a related 77' mass 
modification. This model and the related reduction of the effective intercept parameter of the 
two-pion Bose-Einstein correlation function was proposed first in ref. [16) and detailed recently 
in Refs. Q31H1HS]. 

In thermal models, the production cross sections of the light mesons are exponentially 
suppressed by the mass. Hence one expects about two orders of magnitude less 77' mesons from 
the freeze-out than pions. This suppression, however, may be moderated as a consequence of 
a possible 77' mass reduction, and the 77' mesons may show up in an enhanced number. The 
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Figure 1: Datascts of A*(mx) from RHIC 
y'spfN — 200 GeV like-sign pion correlation 
measurements from Refs. [TU1 ITT1 IT2"] . 



2 



Vcrtcsi, Csdrgo, Sziklai - In-medium reduction of the rf mass 



number of in-medium rf mesons is calculated with an improved Hagcdorn formula yielding the 
following Tj enhancement factor: 



This formula includes a prefactor with an expansion dynamics dependent exponent a RJ 1 — d/2 
for an expansion in d effective dimensions [20]. As a default value, a — was taken [16j 
and, for the systematic investigations, this parameter was varied between —0.5 < a < 0.5. 
Other model parameters and their investigated ranges are described as follows: T con d in the 
above formula corresponds to the temperature of the medium when the in-medium modified 
if mesons are formed; its default value was taken to be T con d = 177 MeV [H] and varied 
systematically between 140 and 220 MeV. Resonances with different masses were simulated 
with a mass dependent slope parameter T e ff = Tpo + m(uT) 2 , where the default values of 
Tfo = 177 MeV and (ut) = 0.48 [UJ were utilized and systematically varied in the range of 
100 MeV < T FO < 177 MeV and 0.40 < (u T ) < 0.60 . 

Once produced, the rj' is expected to be decoupled from other hadronic matter, since its 
annihilation and scattering cross sections are very small [7J. If the 7/ mass is reduced in the 
medium, the observed 77' spectrum will consist of two components. If the px of the ?/ is 
large enough, it can get on-shell and escape. This will produce a thermal component of the 
spectrum. Energy conservation at mid-rapidity implies m* 2 + Pt r)' 2 = m r)' 2 + Pt,i)' 2 - (In 
the latter equation the quantities marked with an asterisk denote the properties of the in- 
medium 77', while the ones without an asterisk refer to the free 77'.) On the other hand, rj'-s 

with p* T r]l < to*, 2 — rrirji 2 will not be able to leave the hot and dense region through thermal 
fluctuation since they cannot compensate for the missing mass [7j [8] , and thus will be trapped in 
the hot and dense region until it disappears. As the energy density of the medium is dissolved, 
the effect of QCD instantons increases and the trapped 77' mesons regain their free mass and 
appear at low pt ■ 

The low pt enhancement of the 77' meson also affects the spectrum of its decay products, 
first of all the 77, the feed-down to which can be described with an 77' to 77 branching ratio 
of BR(rj' — > 77 + tttt) rj 65.7%. A connection between the 77' enhancement f V ' and the 77 
enhancement f v can be expressed as 



where and N v denote the resonance multiplicities of the 77' and 77 as fixed by the input 
model. 

In our recent works of Rcfs. [T3J [TJJ [T5], we improved on earlier simulations of Ref. (IT)] , 
that considered the trapped 77' mesons to leave the dissolving medium with a negligible ■ 
That earlier approach resulted in a steep hole in the extrapolated intercept parameter A*(tot) 
at a characteristic transverse mass of tot < 250 MeV [Til HH1 H?] . In that simplified scenario 
the only free parameter was the in-medium 77' mass, determining the depth of the observed 
hole. In our recent analysis [13j [TU [15] , the rj'-s from the decaying condensate were given a 
random transverse momentum, following Maxwell-Boltzmann statistics with an inverse slope 
parameter -B -1 , which was necessary to obtain a quality description of the width and the slope 
of the A*(m T ) data of PHENIX and STAR in the to t rj 300 MeV region. Physically, B^ 1 is 




(1) 



1 V / 

fv = 1 + (./V - 1) tt br (v' + ™) 



(2) 
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limited by Tfo, so the trapped rf -s may gain only moderate transverse momenta. Hence, the 
enhancement mostly appears at low pt [6j El [8J just as in the first simulations. However, now 
the slope of "hole" of the A* (tot) curve is determined by and, for certain values of the 

model parameters, the data can be reproduced quantitatively. (The A* values, actually used 
in the presented analysis, and their total errors are discussed in details in Ref. [15]. Here A™ ax 
is the A*(tot) value taken at hit = 0.7 GeV, with the exception of the STAR data, where the 
data point at the highest tot = 0.55 GeV is considered. Note that the tot dependency of the 
A* (tot) measurements in the 0.5-0.7 GeV region is very weak.) 




Figure 2: The transverse mass dependence 
of the relative intercept parameter in the 
PHENIX and STAR dataset is reproduced 
with an in-medium mass modification of the 
rf mesons using two different resonance mod- 
els as input. The same resonance models, but 
without in-medium mass modification, cannot 
explain these datascts. AMPT is known to be 
fairly successful in describing the HBT radii 
within a non-thermal scenario and without rj' 
mass modification [29], however, it is not capa- 
ble of describing the current dataset in a sta- 
tistically acceptable manner. 



Figure 3: Standard deviation contours on the 
(i? -1 , to*,) plain, obtained from A st (TO T )/A™ ax 
of Monte Carlo simulations based on parti- 
cle multiplicities using two different models 
for hadronic resonances. The region between 
the horizontal solid lines indicates the theoret- 
ical range predicted by quark model consider- 
ations, while the dotted horizontal line stands 
for Weinberg's lower limit. 



We have investigated a broad class of models of resonance production, including FRITIOF [22 
and UrQMD [53] models, that both produce resonances without assuming local thermalization. 
Resonance decays, including decay chains, were simulated with JETSET 7.4 [24] . 

The FRITIOF [22] Monte Carlo model, which is based on the superposition of nuclcon- 
nuclcon collisions and the Lund string fragmentation model, cannot describe the behavior seen 
in A» (tot) / A™ aa: even when an arbitrary rf mass modification is considered. On the other hand, 
hadronic cascade based UrQMD [23], as well as the quark coalescence model ALCOR [25] and 
the thermal resonance production models of Refs. [IHl [27] [25], provide a successful fit in a 
certain range of the in-medium rf masses. The main difference between the thermal models 
that we utilized was in those resonance multiplicities that are not yet measured well: ref. |26j 
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predicts a factor of 1.6 more 77-s and a factor of 3 more rj'-a than the models of ref. (27] [28] . 
The relevant resonance fractions of these models are detailed in Table V of ref. [15] . 

The dotted line in Fig. [2] indicates a scenario without an in-medium 77' mass reduction, while 
the dot-dashed and solid lines show the enhancement required to describe the dip in the low 
tut region of A* corresponding to the resonance multiplicities of Refs. EXI ^ respectively. 

Based on extensive Monte-Carlo simulations, x 2 of the fits to the data of Fig.[2]was computed 
as a function of m*, and B^ 1 for each resonance model and each fixed value of model parameters 
of a, T conc [, Tpo and (ut). The best values for the in-medium mass of 77' mesons are in, or 

slightly below, the range <J ^{2rn^ + m 2 .) < m*, < ■\/2m 2 c — to 2 , predicted in ref. [7], while all 

are above the lower limit of m*, > y/im^ given by ref. [30]. The A»(mT)/A™ al simulations for 
the best fits of two characteristic models are compared to the no-mass-drop scenario on Fig. [2] 
while the 1, 2 and 3-a parameter boundaries are indicated in Fig. [3] Models that describe 
both PHENIX and STAR A»(mT)/A™ aa: data in a statistically acceptable manner with the 
assumption of a sufficiently large in-medium 77' mass reduction are all used for the estimation 
of systematics. The key parameters of the best fits are listed in Table [T] 

3 Results 

We have used different input models and setups to map the parameter space for a twofold 
goal: to determine, at least how big 77' in-medium mass reduction is needed to be able to 
describe these datasets, and also to determine, what are the best values of the in-medium mass 
modification of the 77' mesons. Utilizing our indirect method, we have also reconstracted the 
transverse mass dependent spectrum of these 77' mesons. 

3.1 Lower limit on the in-medium rf mass reduction 

We excluded certain regions where a statistically acceptable fit to the data is not achievable, 
thus we can give a lower limit on the 77' mass modification. At the 99.9 % confidence level, 
corresponding to a more than 5-cr effect, at least 200 MeV in-medium decrease of the mass 
of the 77' (958) meson was needed to describe both STAR 0-5 % central and PHENIX 0-30% 
central Au+Au data on A,(jtit) / '\™ ax in y/s^N = 200 GeV Au+Au collisions at RHIC, in the 
considered model class. 

3.2 Best value of the in-medium rf mass reduction 

We have determined the best values and errors of the fitted m*, and B~ l parameters. The 
best simultaneous description of PHENIX [TU] and STAR fTJ] relative intercept parameter data 
is achieved with an 77' mass that is dramatically reduced in the medium created in central 
Au+Au collisions at RHIC from its vacuum value of 958 MeV to 340 +g° + 2 ®° ± 45 MeV. The 
first error here is the statistical one determined by the 1-a boundaries of the fit. The second 
error is from the choice of the resonance model and the parameters (a, T conc i, Tpo and (ut)) 
of the simulation. The third error is the systematics resulting from slightly different PHENIX 
and STAR centrality ranges, particle identification and acceptance cuts. These effects have 
been estimated with Monte-Carlo simulations, detailed in ref. [T5], not to exceed 9.8%, 7% and 
3% respectively. The main source of systematic errors is the choice of the resonance models. 
This is due to the unknown initial 77' multiplicity, hence models like ref. [26] with larger initial 
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rj abundances require smaller in-medium rj mass modification, as compared to the models of 
ref. [28JE7]. 



3.3 Transverse mass spectra of rj and rj mesons 

In addition to the characterization of the in-medium rj mass modification, the transverse mo- 
mentum spectra of the r\ and rj mesons have also been reported in Ref. |13j . 

The reconstructed spectrum of rj and rj in ^/snn = 200 GeV Au+Au collisions, shown for 
selected models in Figs.0]and[5]respectively, feature the characteristic low transverse momentum 
enhancement discussed above. Normalization was carried out with respect to the rj' multiplicity 
of the model described in Ref. [55] . Although PHENIX measured before the rj spectrum in the 
Pt > 2 GeV region [31], as far as we know the spectrum of the rj particles has not been 
determined before in t/snJj = 200 GeV Au+Au collisions at RHIC. Let us note that the 
enhancement of the rj production affects the pt < 1 GeV region only. The higher part of the 
rj spectrum serves as a consistency check when compared to more direct measurements [15] . 




Figure 4: The transverse mass dependent 
spectrum of the rj mesons, obtained using two 
different resonance models as input. The band 
indicates the systematic error, obtained from 
varying the resonance models as discussed in 
the text. 



Figure 5: The transverse mass dependent 
spectrum of the ?/ mesons, obtained using two 
different resonance models as input. The band 
indicates the systematic error, obtained from 
varying the resonance models as discussed in 
the text. 



4 Discussion 

The mass drop effect has been inspected from several aspects including the dependence on 
the properties of the colliding systems. Note however, that the validity of our analysis relies 
on the correctness of published STAR and PHENIX data. These data, however, were not 
measured with the definite purpose to serve as a base for the search for the in-medium rj 
mass modification, where particular attention has to be payed to the momentum dependence 
of the particle identification purity and efficiency, especially at low px regime. More detailed 
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Resonance 


Best fit to combined STAR+PHENIX dataset 




5-(T limit on 


model 


m*, (MeV) 


S" 1 (MeV) 


X 2 /NDF 


CL % 


f* 


fn 


m*, (MeV) 


ALCOR [25] 


490 ™ 


42 


20.2/11 


4.29 


43.4 


5.25 


< 700 


Kaneta [55] 


530 


55 


22.8/11 


4.12 


25.6 


3.48 


< 730 


Letessier [2"T] 


3 4 q+50 


86 


18.9/11 


6.35 


67.6 


4.75 


< 570 


Stachel [28] 


340 +5 ° 

° -60 


86 


18.9/11 


6.38 


67.6 


4.97 


< 570 


UrQMD [23] 


400 +50 


86 


18.9/11 


6.14 


45.0 


7.49 


< 660 



Table 1: Best fits of m* , and B 1 for different input resonance multiplicity models, followed by 
X 2 1 NDF and the corresponding confidence level (CL) and the integrated rf and r\ enhancement 
factors frf' and /,,. The 5-cr limits of maximum in- medium masses are also shown. Errors on 
m*, values represent l-cr boundaries, while the 5-cr limits include systematic errors too. 



dedicated A H ,(mx) measurements together with additional analysis of the dilcpton and photon 
decay channels of rf could help consolidate the findings reported here. 



4.1 System size, centrality and energy dependence 
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Detailed analysis of the STAR and PHENIX 
A*(m T )/A™ ax dataset recorded at 7.7, 9.2, 
11.5, 39 and 62.4 GeV during 2010 has just 
been started [32j , marking the beginning of the 
RHIC energy scan program. 

Detailed data available from the NA44 col- 
laboration at sjv n = 19.4 GeV [TI5] as well as 
from the STAR collaboration at different cen- 
tralities within the 0%-80% range at ^snn = 
62.4 and 200 GeV Cu+Cu and Au+Au colli- 
sions [18] is shown on Fig. [6] The NA44 data 
at -y/ijvw = 19.4 GeV does not feature an rf 
mass drop effect. However a positive sign of 
the rj 1 mass modification is apparent in each 
case of the STAR datascts, indicating that the 
mass modification effect is nearly at maximum 
in y/SNN = 200 GeV Au+Au collisions and re- 
duces with decreasing centrality, colliding en- 
ergy and system size. We have estimated the 
magnitude of the system size and energy de- 
pendence between 62.4 GeV Cu+Cu and 200 

GeV Au+Au collisions to be not larger than 15%, which is substantially less than the dominant 
systematic error coming from the choice of the resonance model. 



A.(k T )/X maJ( comparison 
■ NA44 Pb+Pb 
i STAR Au+Au 62GeV 0-5% 

* STAR Au+Au 200GeV 0-5% 
STAR Cu+Cu 62GeV 0-10% 

* STAR Cu+Cu 200GeV 0-10% 



100 150 200 250 



300 350 400 

k T [MeV] 



450 500 550 600 



Figure 6: The energy and system size depen- 
dence of the relative intercept parameter in the 
NA44 S+Pb and STAR Cu+Cu and Au+Au 
datasets. 
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4.2 Channels other than BEC 



A promising channel of observation is the dilcptonic decay 77' —> £ + £~-f, because a low-px 
77' enhancement would give extra lepton pairs to the low invariant mass region. The paper 
of Kapusta, Kharzeev and McLerran on the return of the prodigal Goldstone boson [7J was 
in fact motivated by the dilepton enhancement seen in CERN SPS 200 AGeV/c energies in 
S+Pb reactions. Recent interpretations of CERES [33J and NA60 data [3J] indicate that the 
approach to a chiral symmetry restored state could proceed through resonance broadening 
and eventually subsequent melting, rather than by dropping masses or mass dependency or 
mass degeneracy between chiral partners [35]. Recent PHENIX findings also show a definite 
excess in the m e + e - < 1 GeV dielectron invariant mass region in -^/SaT/v = 200 GeV Au+Au 
collisions [36) . Unlike at lower beam energies, in this case the contribution from a hot hadronic 
phase without mass shifts seems to be insufficient to account for the enhancement seen in the 
data [37] . Estimations using the enhancement factors in Table Q] indicate that the observed 
in-medium rf mass drop is indeed a promising candidate to explain this dilepton excess. 

PHENIX recently reported a two-component transverse momentum spectrum in dilepton 
channel direct photon measurements [36] . which provides an additional testing possibility to 
constrain the two component structure of the 77' spectra reported here. 

5 Summary 

Our report presents a statistically significant, indirect observation of an in-medium mass mod- 
ification of the 77' mesons in ^/snn = 200 GeV Au+Au collisions at RHIC. These results were 
recently published in Refs. |13l H4l [T5]. A similar search for in- medium 77' mass modification 
provided negative result in S+Pb reactions at CERN SPS energies [IB]. More detailed studies 
of the excitation function, the centrality and system size dependence of the relative intercept 
\*{mT)/\™ ax could provide important additional information about the onset and saturation 
of the partial J7a(1) symmetry restoration in hot and dense hadronic matter. Studies of the 
low-mass dilepton spectrum and measurements of other decay channels of the 77' meson may 
shed more light on the reported magnitude of the low pt 77' enhancement and the related UaO-) 
symmetry restoration in high energy heavy ion collisions. 
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